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CNJ | Abstract 

In this article, we analyze a two-level finite element method for the two dimensional time- 
t ^~ ■ dependent incompressible Navier-Stokes equations with non-smooth initial data. It involves solv- 

ing the non-linear Navier-Stokes problem on a coarse grid of size H and solving a Stokes problem 
on a fine grid of size h, h « H. This method gives optimal convergence for velocity in ii^-norm 
and for pressure in L 2 -norm. The analysis takes in to account the loss of regularity of the solution 
at t = of the Navier-Stokes equations. 
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1 Introduction 
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In this article, we consider a two-level semi-discrete finite element approximation to the two dimen- 
sional time-dependent incompressible Navier-Stokes equations 

(1.1) y + U'Vu - i/Au + Vp = f(x, t), x G n, t > 
I with incompressibility condition 

(1.2) V • u = 0, x G O, t > 0, 

•i-H ■ 

and initial and boundary conditions 

(1.3) u(x,0) = uo in fi, u = 0, on dQ, t>0. 

Here, is a bounded domain in M. 2 with boundary dQ and v > is the viscosity, u and p are the 
velocity field and pressure, respectively. And f is a given force field. 

Two-level or two-grid methods are well-established and efficient methods for solving non-linear 
partial differential equations. Due to high computational cost of solving a non-linear problem on 
a fine grid, we solve the original problem on a coarse grid and update the solution by solving a 
linearized problem on a fine grid. In other words, in the first step, we discretize the non-linear PDE 
on a coarse mesh, of mesh-size H and compute an approximate solution, say, u#. Then, in the 
second step, we formulate a linearized problem, out of the original one, using u# and discretize it 
on a fine mesh, of mesh-size h, h << H, thereby, compute an approximate solution, say, u h . With 
appropriate h,H, we obtain same order of convergence of the error u — u h , as that of the error 
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obtained by semi-discrete finite element Galerkin approximation on the find grid; but with far less 
computational cost, since, instead of solving a large non-linear system, we solve a small non-linear 
system and a large linear system. 

In this article, we study the following two-level finite element approximation for the problem 
(|l.l|) - (|1.3p : First, we compute a semi-discrete Galerkin finite element approximations {vlhiPh)i ° ver 
a coarse mesh of mesh-size H. Then, we use the approximation uh to compute a semi-discrete 
Galerkin finite element approximations (u h ,p h ) of the following Stokes problem: 

(1.4) U(-i/Au + u/ I -Vu fl + Vp=f, V-u = 0, in ft 

over a fine mesh of mesh-size h « H. 

The above algorithm is nothing new and in fact, this and similar algorithms have been studied 
on numerous occasions. But to the best of our knowledge, no study has been done for non-smooth 
initial data. Our main aim in this work is to do error analysis of the above mentioned two-level 
method, under non-smooth initial data. 

Two-grid method was first introduced by Xu [23} [24"] for semi-linear elliptic problems and by 
Layton et.al [HI H5j HB] for steady Navier-Stokes equations. It was carried out for time dependent 
Navier-Stokes by Girault and Raviart [7J for semi-discrete case. The method may vary depending 
on the algorithm, by formulating different linearize problem, to solve in the second step; like, in the 
case of Navier-Stokes, one can chose a Stokes problem or an Oseen problem or a Newton step to 
solve on the fine mesh. Several works in this direction, involving both semi-discrete and fully discrete 
analysis, can be found in [H [21 El QUI OH EH [EH QH [20] and references therein. 

The two-level methods are similar to non-linear Galerkin methods, post-processed and dynamical 
post-processed methods, in the sense that, all these method try to control the computational cost 
or efficiency by controlling the non-linear term, implementing similar ideas: solve the non-linear 
problem on a coarse grid and use that solution to solve a linearized problem on a fine grid. A 
discussion on this can be found out in [5]. 

But in all these methods, including two-level, very few articles can be found to carry out the 
analysis, taking into account the loss of regularity of the solution of the Navier-Stokes equations at 
the initial time. Recently, two articles [H [6] have taken this into account. For example, both these 
articles assume no more than second-order spatial derivative of the velocity bounded in L 2 , up to 
initial time t = 0. Under realistic assumptions on the given data, the solution of the Navier-Stokes 
equations suffers from lack of regularity at t = 0. For us to assume regularity at t = 0, we need that 
the data satisfy some non-local compatibility conditions at t = (see [E]). These conditions are 
very hard to verify even for simple models and hardly appear in physical context. We have avoided 
these conditions in our analysis and more. We have worked here with non-smooth initial data, that 
is, initial velocity uo belong to Hq. Now, even the second-order spatial derivative of the velocity is 
not bounded in L 2 , up to initial time t = 0. Lemma 12.11 tells us that ||u(i)||2 ~ 0(i -1 / 2 ). So, our 
assumption on the initial data forces us to have only the first-order spatial derivative of the velocity 
bounded in L 2 at t = 0. This proves to be a bottle-neck in our error analysis. 

Our work here closely resembles the work of [9j and we present here similar results that have 
been obtained in [9], but with weaker condition on uo. In that sense, this work is devoted to the 
technicalities required to deal with the non-smooth initial data. 

Apart from the error analysis of this two- level method, we have also discussed the error analysis 
of finite element Galerkin method for (jl.ip - (|1.3j) . Although, these later results are discussed in 
[12] . the results here vary slightly, due to non-smooth initial data. Also, the pressure estimate has 
been improved for conforming finite element, in the sense that, it now reads as ~ 0(t _1 / 2 ) rather 
than ~ 0(t _1 ). For smooth initial data, it can be shown to have no singularity at t = (using 
similar argument presented here, e.g., using Lemma |4. II for pressure error), rather that the estimate 
~ 0{t~ l l 2 ) in [12] , Similar improvement is also realized in the two-level method. Theorem 5.3 in 
[9] states that the error in velocity is ~ 0{H 2 ) and the error in pressure is ~ 0{t- 1 / 2 H 2 ), when the 
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initial data is smooth. Note that whereas the velocity error remains bounded at t = 0, it is not so for 
the pressure error. In our case, we have shown that velocity and pressure errors are ~ 0(i _1 .ff 2 ) (see 
Lemmas 16.31 and 16. 4|) . when initial data is non-smooth. By following similar procedure, for smooth 
initial data, we can show that the errors, in fact, are ~ 0(H 2 ). This improves the presure error 
presented in [9]. 

Finally, we observe that it is of practical importance to discretize in time and to implement the 
fully discretize scheme to verify the usefulness of the present method. In fact, computational results 
are shown in [11] to support the fact that this two-level method is better that the standard Galerkin 
finite element method. But the analysis presented there, for a fully discrete two-level method is for 
smooth initial data, that is, the initial velocity uo is in Hq n H 2 . We will deal with non-smooth 
initial data for a fully discrete two-level method somewhere else. 

In the following sections, we will assume C and K to be generic positive constants, where K 
generally depends on the given data, that is, uo,f. And for the sake of convenience, we would write 
v(t) as v when there arises no confusion. 

The article is organized as follows. In section 2, we introduce some notations and preliminaries, 
and present a priori estimates and regularity results for the solution of (jl.ip - (jl.3j) . when the initial 
data is non-smooth. Section 3 deals with finite element Galerkin approximations, whereas Section 4 
deals with the error analysis of the Galerkin approximation. The two-level method is briefly discussed 
in Section 5 and its error analysis is carried out in Section 6. 



2 Preliminaries 

For our subsequent use, we denote by bold face letters the Revalued function space such as 

Ho = [Hq(Q)] 2 , L 2 = [LW. 
Note that Hj is equipped with a norm 

2 2 

l|Vv|| = {Y j {d j v i ,d j v i )) 1/2 = (^(V^,V^)) 1/2 . 

i,j=l i=l 

Further, we introduce some divergence free function spaces: 

Ji = {</> G Hq : V • = 0} 
J = {(p £ L 2 : V • 4> = in Q,(f> ■ u\qq = holds weakly}, 

where n is the outward normal to the boundary d£l and 4> ' n |c9f2 = should be understood in the 
sense of trace in H" 1 / 2 ^), see [2T]. For any Banach space X, let L p (0,T; X) denote the space of 
measurable X -valued functions <j> on (0, T) such that 



I \\4>(t)\\ p x dt < oo ifl<p<oo, 
J o 



and for p = oo 

ess sup \\4>(t) \\x < oo if p = oo. 

0<t<T 

When there arises no confusion, we simply denote these spaces by L P (X). Through out this paper, 
we make the following assumptions: 

(Al). For g € L 2 , let the unique pair of solutions {v € J\,q € L 2 /M.} for the steady state Stokes 
problem 

-Av + Vq = g, 
V • v = in Q, v\ 8 n = 0, 
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satisfy the following regularity result 

IMb + IMIifi/K < c||g||. 

(A2). The initial velocity uo and the external force f satisfy for positive constant Mq, and for T 
with < T < 00 

u G Ji, f,f t ,f tt G L°°(0,T;L 2 ) with ||u ||i < M , sup {||f||, ||f t ||, ||f tt ||} < M . 

0<t<T 

With P : L 2 — J as orthogonal projection and A = P(— A) : Ji n H 2 as Stokes operator, we first 
note that (Al) implies 

||v|| 2 < C||Av||, veJiPH 2 . 

And (see [12, (2.4)]) 

(2.1) ||v|| < A^ 1/2 ||v||i < Ai 1 ||v|| 2 , v G Ji n H 2 , 
where Ai > to be the least eigenvalue of the Stokes operator. 

Before going into the details, let us introduce the weak formulation of (|l.lj) -( lL3]l . Find a pair of 
functions {u(t),p(t)}, t > 0, such that 

(2.2) (u t) 0) + I /(Vu,VtfO + (u-Vu,0) = (p,V-0) + (f,0) V0 G Hj, 

(V-u, X ) = V X GL 2 . 

Equivalently, find u(t) G Ji, t > such that 

(2.3) (u tj cf>) + i/(Vu, V0) + (u • Vu, 0) = (f, 0), V0 G Ji. 

For existence and uniqueness and the regularity of the solution of the problem (12.2p or (12. 3p . we 
refer to [12]. For the sake of completeness, we present below the a priori estimates and higher-order 
estimates, which will be used in our error analysis. These vary slightly with the results presented in 
[12j due to non-smooth initial data. And hence, we sketch a proof here. 



Lemma 2.1. Assume that (Al) and (A2) hold and let < a < i>\i. Then, for some constant 
K > 0, which depends only on the given data, the weak solution pair (u,p) of il.l\) - [L3\) satisfies 
the following estimates: 



(2.4) 



\\u(t) || 2 + e~ 2at f* e 2as ||u(a) \\jds < K, 
Jo 

(2.5) \\u(t)\\i + e- 2at r e 2QS {||u( S )|| 2 + ||u s ( S )|| 2 } ds<K, 

Jo 

(2.6) r*(t){||u,(t)|| 2 + ||u(t)|| 2 + \\p(t)\\l} + e~ 2at f a(s)\\u s ( S )\\ 2 ds < K, 

Jo 

(2.7) (r*(t)) 2 \\u t (t)\\ 2 + e- 2at [\ 1 (s){\\u s (s)\\ 2 + \\uMf + \\Ps(s)\\l} ds < K, 

Jo 

(2-8) (r*(t)f/ 2 {\\u tt (t)\\ + ||u t (t)|| 2 + ||p t (t)|| 2 } < K, 

where r*(t) = min{l,t}, <r(t) = T*(t)e 2crf and (7i(i) = (t* (t)) 2 e 2at . 
Proof. We choose = e 2at u = e at u in (I2.3P and use (12. lj) to find 

(2.9) --llup+^-^IH^—llfllllull, 
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Use kickback argument, then integrate and multiply by e 2at to obtain 

u| ' 2 + ( - - B e ' 2 " II l|i(s)ll?<is 5 e " 2 " ,||u » 112 + »S ||f|l » = : * 

where ||f ||oo = ||f ||l°°(l 2 )- 

Now choose (f> = u in (12. 3p and integrate from t to t + to for some fix to > 0. 

(2.10) ||u(t + t )|| 2 + [ t+t ° \Hs)\\lds < ||u(t)|| 2 + ^° ||f \\l < K(t ). 

Jt v 

Next, we choose 4> = Au in J23D to find (see [12 (2.11)]) 

(2.11) ^||u|| 2 + z,||Au|| 2 < C||f|| 2 + C||u|| 2 ||u||f. 

Keeping in mind (|2.10p . we apply uniform Gronwall lemma (see \12\ Lemma 1.1, pp. 91]) to observe 
that for fixed to > 

||u(t + t )||i <K, t > 0. 

Multiply (12. lip by e 2at , integrate with respect to time, use the above estimate to conclude the first 
part of (|2.5p . Choose = e 2at u t for the rest of (|2.5p . 

For next estimate, we differentiate (|2.3p with respect to time and put <fi = cr(t)u t to get 
^{a(t)\\u t \\ 2 } + 2wr(t)||u i || 2 = a t (t)\\u t \\ 2 - a(t)(u t • Vu, u t ) + <r(t)(f t , u t ). 
Observe that 

(u t • Vu,u f ) < ||u t ||L4(Q)||u||i < C||u t ||||ut||i||u||i. 

And therefore, we have 

^{a(t)||ui|| 2 } + ^(t)||u t || 2 < Ce 2at \\u t \\ 2 (l + ||u|| 2 ) + a(t)||f t || 2 . 

Integrate with respect to time and use (|2.5p to conclude part of (|2.6p . The rest of it can be proved 
by using the equations (|2.3p and (|2.2p . 

Put 4> = (Ji(t)Aut and = o"i(t)u tt , respectively, after differentiating (|2.3p and proceed in similar 
fashion as above to find (<Ji(i) = (r*(t)) 2 e 2Qi ) 

(2.12) (r*(t)) 2 ||u t (t)|| 2 + e- 2at /" <n(s){ ||u s ( S ) || 2 + ||u ss ( S ) || 2 } ds < K. 

Jo 

Differentiate the equation (|2.2p and use f|2. 12j) to obtain the pressure estimate of (|2.7p . 

Finally, we differentiate (|2.3p twice, with respect to time. First, put <f> = o"2(t)u#, 02 (t) = (r*(i)) 3 e 2a * 

and use (|2.7p to obtain 



(r*(t)) 3 ||u ti (t)|| 2 + e- 2at /V 2 (s)||u ss ( S )|| 2 ck< 

•/ o 



if. 



And then use the double differentiated equation and the equation (|2.2|) (after double differentiation) 
with the above obtained estimate to conclude (j2.8|) and this completes the rest of the proof. □ 

We present below Gronwall's Lemma, which will be used subsequently. 
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Lemma 2.2 (Gronwall's Lemma). Let g,h,y be three locally integrable non-negative functions on 
the time interval [0, oo) such that for all t > 

y(t) + G(t) < C + [ h(s) ds+ I g(s)y(s) ds, 
Jo Jo 

where G(t) is a non-negative function on [0, oo) and C > is a constant. Then, 

y(t) + G{t)< (c+ J h{s) ds^exp^J^ g(s) ds^. 

3 Classical Galerkin Method 

From now on, we denote h with < h < 1 to be a real positive discretization parameter tending 
to zero. Let and L^, < h < 1 be two family of finite dimensional subspaces of Hq and 
L 2 /M, respectively, approximating velocity vector and the pressure. Assume that the following 
approximation properties are satisfied for the spaces and L^: 

(Bl) For each w G Hq Pi H 2 and q G H 1 /R there exist approximations ihW G and jhq G L^ such 
that 

||w - i h w\\ + /i||V(w - ihw)|| < Ch 2 \\w\\ 2 , \\q-jh<l\\ < C7*IM|i- 
Further, suppose that the following inverse hypothesis holds for w/j G H^: 

(3.1) ||Vw fc || < Ch^WwhW. 
To define the Galerkin approximations, we set for v, w, cb G Hq, 

o(v,0) = (VV.V0) 

and | 

6(v,w,0) = -(v • Vw,<£) - -(v • V0, w). 

Note that the operator &(•,-,■) preserves the antisymmetric property of the original nonlinear term, 
that is, 

H v h, w fc , w ft ) = Vv^, w ft G H h . 

The discrete analogue of the weak formulation (|2.2p now reads as: Find u/j(i) € H^ and Ph(t) G L/, 
such that u/ l (0) = uq^ and for t > 

{u ht ,<t> h ) + va(xi h ,(t> h ) + 6(u ft , u/j, - (p ft , V • h ) = (f , <t> h ), 

(3.2) (V-Ufc.x/,) = 0, 

for fe G H/j, Xfc ^ Here u h G H/j is a suitable approximation of uo G Ji. 

In order to consider a discrete space analogous to Ji, we impose the discrete incompressibility 
condition on H/j and call it as J^. Thus, we define J^, as 

Jft = {vh G H h : (xh, V • v A ) = Vxh G L/J. 

Note that is not a subspace of Ji. With as above, we now introduce an equivalent Galerkin 
formulation as: Find Uh(t) G Jh such that u/ l (0) = uoh and for t > 

(3.3) (u ht ,(j) h ) + ua(u h , (f> h ) = -b(u h ,u h ,(f> h ) + (f,</> ft ) V0 h G J ft . 
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Since Jh is finite dimensional, the problem (|3.3p leads to a system of nonlinear ordinary differential 
equations. For global existence of a unique solution of (|3.3j) (or unique solution pair of (|3.2p ). we 
again refer to [TZ] . 

For continuous dependence of the discrete pressure ph(t) G Lh on the discrete velocity Uh{t) G Jh, 
we assume the following discrete inf-sup (LBB) condition on the finite dimensional spaces H/j and 
L h : 

(B2') For every G Lh, there exists a non-trivial function <p h G H/j and a positive constant i^O; 
independent of h, such that 

|(?fc,V-0 h )| > C\\V<P h \\\\q h \\ L 2 /M . 

Moreover, we also assume that the following approximation property holds true for Jh- 
(B2) For every w G Ji H H 2 , there exists an approximation r/,w G Jh such that 

||w — r/jw|| + /i|jV(w — r^w)!! < C/i 2 ||w||2. 

The I? projection Ph : L 2 !->■ satisfies the following properties (see [12]): for G Jh, 

(3.4) 110-P^II + ^HVP^II < C^||V0||, 
and for G Ji l~l H 2 , 

(3.5) ||0 - P h <t>\\ + h\\V{<j> - P h <t>)\\ < a^ 2 ||A0||. 

We now define the discrete operator : i— >• through the bilinear form a(-, •) as 

(3.6) a(v h , 4>h) = (-AftVft, 0) Vvft, <p h G H^. 

Set the discrete analogue of the Stokes operator A = P(— A) as A^ = Ph(-Ah). Note that the Ah 
restricted to Jh is invertible and we denote its inverse by (Ah)- 1 ; for details, see [13 [p3]. Following 
[13] . we define discrete Sobolev norm as: 

K|| r = \\(-A h ) r/2 Vh\\, for v h G Jh, T G M. 

Using Sobolev imbedding and Sobolev inequality, it is easy to prove the following Lemma (similar 
to m (3-4)]). 

Lemma 3.1. Suppose conditions (A.1), (Bl) and (B2) are satisfied. Then there exists a positive 
constant K such that for v, w, G H^, the following holds: 

|v|| 1 /2|| V vl[ 1 /2|| Vw |]l/2|| AftW ||l/2||^ ) j| 5 

IvIl^llA^vll^UVwIII^II, 

|v|| 1 /2||Vv|| 1 /2||Vw||||0|| 1 /2|jV0|| 1 /2 ) 

IvIIHVwIIIMlVailA^llVa, 
|v||||Vw|| 1 /2||A ft w||V2||0||V2||v0||V2 

Examples of subspaces and Lh satisfying assumptions (Bl) and (B2') are abundant in literature, 
for example, see [31 HI [8]. 

We present below a couple of lemmas, one dealing with a priori and regularity estimates of and 
the other, with higher-order regularity results. The proof is similar to that of Lemma 12. 11 

Lemma 3.2. Under the assumptions of Lemmas \2. 1\ and \3. 1\ the semi-discrete Galerkin approxi- 
mation u/i of the velocity u satisfies, for t > 0, 

(3.8) |M*)|| 2 + e- 2a * f e 2as \\u h (s)\\ 2 ds < K, 

Jo 

(3-9) \\u h (t)\\l + e- 2at f e 2as {\\u h (s)\\l+\\u h , s (s)\\ 2 } ds < K, 

Jo 

(3-10) (T*(*)) 1/2 {K(t)|| 2 + \\p\\m/n} < K, 



(3.7) |(v-Vw,0)| <K< 
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where r*(t) = min{l,t} and K depends only on the given data. In particular, K is independent of h 
and t. 

Lemma 3.3. Under the assumptions of Lemma \3.S\ the semi-discrete Galerkin approximation Uh 
of the velocity u satisfies, for t > and for r G {0, 1}, i G {1, 2}, r + i < 2, 



(3.11) (T*(t)) r+2 *- l \\Dlu h (t)\\ 2 r + e- 2at !\r*(s)Y^ i - 1 e 2a ^\D\vi h {s)t +1 ds < K, 

Jo 

where D\ = J^-. And 

(3-12) e~ 2at f\r*(s)fe 2as {\\u Kss (s)\\ 2 + llp^OOIIflV*} ds < K, 

(3-13) (r*(t)) 3 / 2 {K t || 2 + ||p M (t)||i} <K. 



Here K depends only on the given data. In particular, K is independent of h and t. 

4 Error Analysis: Galerkin Method 

In this section, we briefly present a convergence analysis for Galerkin approximation. The analysis for 
smooth initial data, that is, uo G Ji PlH 2 can be found in [12]. The analysis below is for non-smooth 
initial data, that is, uo G Ji. The proofs follow similar lines as those for smooth data (see [12]) apart 
from a few modifications. We will simply try to highlight these modifications in our proofs. 

Theorem 4.1. Let Q, be a convex polygon and let the conditions (A.\)-(A2) and (B1)-(B2) be 
satisfied. Further, let the discrete initial velocity Uq^ € Jh with uq^ = -f\uo, where uq € Ji. Then, 
there exists a positive constant C, that depends only on the given data and the domain £1, such that 
for < T < oo with t G (0, T) 

||(u - u h )(t)|| + fc||V(u - 0(t)|| < Ke K H- l l 2 h 2 . 

Proof. Denoting the Galerkin approximation error as E = u — u^,, we split the error in two parts. 

E = (u - \ h ) + (v h - u h ) =:£ + rj, 

where satisfies the linearized equation (5.3) from [12]. We have (see |12^ Lemma 5.1]; same proof 
will go through) 

(4.1) e~ 2at f ' e 2as U(s)\\ 2 ds<Kh 4 . 

Jo 

For L°°(L 2 ) estimate of £, we now split £ as follows: 

£ = (u - v A ) = (u - S h u) + (S h u - v h ) =:C + 0, 
where Sh is given by (4.52), |12] , Lemma 4.7, [12] tells us that 



(4-2) HCII + h\\C\\i < Ch 2 {(l + ||u||i)||u|| 2 + b||i} 

(4-3) HCtll + h\\Ct\\i < Ch 2 {(l + ||u||i)||u t || 2 + ||p t ||i} 

In order to complete the estimate for £, we only need to estimate 9. The equation in reads as 
(4.4) (6 t ,(j> h ) + va(0,(f> h ) = -(C t ,<t> h ), cf> h £Jh. 
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For a(t) = r*(t)e 2at , we put cf) h = o{t)0 in flHD to find 

j t {v{t)\\e\\ 2 } - a t {t)\\ef + 2ua(t)\\e\\l = -2o(t)(Ct, o). 
Using Cauchy-Schwarz and Young's inequality, we obtain 

±{a(t)\\ef}+2ua(t)\\e\\l < Ce 2at \\0\\ 2 + ai(t)\\Ct\\ 2 

<Ce 2Qt {||e|| 2 + ||CH 2 } + ^(t)IIC t || 2 . 
Integrate with respect to time and use (I4.ip . (|4.2j) and (|4.3p . 



a{t)\\e(t)f + v I a(s)\\e(s)\\lds<Ke 2od h 4 + Ch^J e 2as { \\u(s)g + \\p(s)\\ 2 } ds 

+ Ch* rai( S ){||u s ( S )||| + ||p s ( S )|| 2 } 



Use ([231) and §F7§ and then multiply by e~ 2at to get 

r*(t)\\e(t)\\ 2 + /%)||0( S )|| 2 d S < Kh\ 

Jo 

Now an use of triangle inequality along with inverse hypothesis (|3.ip results in 

(4.5) \W)\\+hM(t)\\l<Kt- 1/2 h 2 . 
It now remains to estimate 77. The equation in 77 is 

(4.6) (»7t,0h) + fa(»7,#h) = Ai,h(# fc ), A €Jfc, 
where 

(4.7) ^x,h{4>h) = H u h,^h, 4>h) ~ K u > u > 0ft) = E, h ) - 6(E, u, <£ A ). 

Due to non-smooth initial data, we need an intermediate estimate, before we proceed for L°°(L 2 ) 
estimate of 77. First we need L 2 (L 2 ) estimate of 77. For that, choose (f)^ = e 2at (— A^)^ 1 r] in (|4.6f) to 
obtain 

(4-8) ^IMI-i " a II^H-i + ^II^H 2 = ^^((-Afc)- 1 !!). 

Here, 77 = e at r\. Using (|4.7p and the definition of £>(•, •, •), we observe that 

A^d-Ah)-^) = - b(u h ,£ + rj, (-Ah)- 1 ^ - b(£ + 77, u, (-Afc)" 1 ^) 

= - i((uft • V)(€ + 77), (-A,)-^) - ^((u fc • V)(-A ?i )- 1 r7, £ + 77) 
(4-9) - + 77) • V)u, (-A*)" 1 ,,) - \m + n) ■ V)(-A,)- 1 r7, u). 

Use Lemma 13.11 to arrive at the following: 

(4.10) -i(K-v)(-A,rV£ + '7) < crKll^Kil^lltillVfhll^di^H + ||^||) 

(4.11) ~((« + n)-v)u,(-A A )- 1 n) < ^ll + NIDIIulliNII^NII 172 

(4-12) + 77) • V)(-A,)- 1 r7, u) < C{M\\ + llr/IDhll^hir^llull^Hull 1 / 2 . 



For the first term on the right-hand side of (|4.9p . we have 

JL 

dx 

2 



(with the notations A = -Jr- and v = (t>i,t>2)) 



((u h -V)(£+ V ),(-A h )- 1 r ] )= V / UM A{(e + ^}((-A /l )- 1 T7) J 

U=i Ja 

2 2 
= -J2 ! A^(e + ^((-A,)- 1 r 7 ) J rfx- £ / UM A{((-A ft )- 1 r 7 ),}(^ + ^ 

i,j=l i,j=l Jn 
= -((V • u ft )(£ + 77), (-A h ) _1 ?7) - ((uh • V)(-A h )- 1 r 7 , (£ + 77)) 
(4.13) < C\\uUU\\ + hlDllnll^llnir 72 + C||u,|| 1 /2|| u ,||V2|| T? ||V2|| r7 ||i/ 2( ||^|| + ii^Hj 

in view of Lemma 13 .11 Incorporate (|4.10p - (|4,13|) in (|4.9p and use Lemma |2. II to find 

A 1 , h ((-A,)- 1 r 7 ) < CU\\\\v\\ y iH\ 1/2 {\\u h \\i + llullx} + C 1 1 »7 II ^? 1 1^7 1 1 3/2 1 1 " 1 1 1 
(4-14) <^H 2 + ^(l|r7ll-i + llel| 2 )- 

Hence, with appropriate e, we obtain from (|4.8p 

|ll^l| 2 -i+^l| 2 = ^(ll^l| 2 -i + lltl| 2 )- 
Integrate and use (|4.1|) . Apply Gronwall's lemma to coclude that 

(4.15) \\v(t)\\-i + e- 2at f e 2as \\-n{s)\\ 2 ds < Ke Kt h A . 

Jo 

Now choose (f> h = cr(t)r) in (|4.6p to obtain 
d 



As in (14.91). we have 



;ft {^WhH 2 }+2^)l|r?|| 2 = ^(t)hl| 2 + 2a(t)A 1 ,,(r ? ). 



M,h(v) = ~b( u h, £ , »7) - 6(€ + *7, u, 77) 
(4-16) = - i((u A • V)£, V) ~ ~(K • V)77, " + »7) • V)u, 77) - + 77) ■ V)t7, u). 



Using Lemma 13. 11 we find 

- ^(((€ + *?) • v)u, V ) - + 77) • v)r/, u) 

• V)u, 77) - • V)t7, u) - ~((tj ■ V)u, 77) - |((t| ■ V)t7, u) 
< Cll^lllllull^llull^ll^ll^ll^ll^ + ll^ll.llulli/^^llV^ 

(4.17) +tf|MIIMIilMli + chu^i^jif || u ||i/2|| u ||V2. 

(4.18) - i((ufc • V)r,,£) < ^Kll^llufcll^llnllil^ll. 

And as for the first term on the right-hand side of (|4.16p . we observe, as in (|4,13D . 

~((u h ■ V)£, 77) = -i((V • u h )£, 77) - i((u fc • V)77, 
(4-19) <C||u,||i /2 ||u /l ||^ 2 ||^||||7,|| 1 / 2 ||r 7 || 1 /2 + C||u,|| 1 /2|| u ,||V 2 
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Incorporating (|4.17j) - ()4.19|) in (|4.16p and using Lemmas 12.11 and 13.21 we obtain 

Ai )h (T/) < kuwwvUt*®)- 1 / 4 + KiMWnh + ||r7|| 1/2 h||f 2 ) 

(4.20) < K||^||||r 7 ||i(r*(t))- 1 / 4 + eh|| 2 + iilT7|| 2 , 
for some e > 0. Therefore, 

2a(t)A lth (rj) < 4^)||t7|| 2 + Ka(t)\\rj\\ 2 + Ke 2at U\\ 2 - 
With appropriate e, we now obtain 

f t {amv\\ 2 } + Mmv\\i<K{Uf + \\f,f}. 

Integrate and use (|4.ip and ()4.15j) to find that 

(4.21) r*(t)h(t)|| 2 + e~ 2at f ||r/(s) || < Ke Kt h\ 

Jo 

With inverse hypothesis (|3.ip . we conclude that 

(4.22) ||^)|| + h\\ V (t)\\ 1 <Ke Kt t- 1 l 2 h 2 . 

This along with (14.51) completes the proof. □ 



Theorem 4.2. Under the assumptions of Theorem \4-l\ and with an additional assumption o/(B2 ), 
we have 

\\{p-p h ){t)\\<Ke Kt t- l l 2 h. 
Proof. Following Lemma 6.1, |12| . and using (Bl), (I2.6P and (B2') along with the observation that 

MA<f>h) <^l|E||i||^!|i 

we obtain 

(4.23) W(p-Ph)(t)\\ L yR < Kt- l ' 2 h + i^||E(i)||i + C\\E t {t)\U th , 

where 

< Et,0 h > 



|Et||-i,/i = sup 



(#0 h eH h \\<t>h\\i 
Taking supremum over a bigger set, we find that 

(4-24) PiU-v* < HEtll-! = sup < ^ t 'f > . 

This is possible due to conforming finite elements. We now have from (|4.23p . using Theorem 14. 1\ 

\\(P-Ph)(t)\\l?/R < Ke Kt t- l l 2 h + C\\E t (t)\\-i- 
We complete the proof by proving the following Lemma. □ 
Lemma 4.1. The error E = u — due to Galerkin approximation satisfies, for t > 

||Et(t)||_i < Ke Kt r x l 2 h. 
11 



Proof. The equation in E reads as 

(4.25) (E t ,<f> h ) + ua(E,<i> h ) = A 1>fe (0J + (p, V • <f> h €J h . 
Now, for any <fi € Hq, we use (|4.25p to find 

(E t , cb) = (Et, P h 0) + (Et, - P^) 

(4.26) = -i/a(E, P ft 0) + A lth (P h cf>) + (p, V ■ P^) + (E*, - 

Using discrete incompressibility condition, ilf 1 -stability of P^ and (Bl), we obtain 
(4.27) 



f -^a(E,P h 0) < CHEllill^lli 

(p,V-P h cj>) = (p - j hP ,V ■ P h cj>) < C%||i||. 
Ai >h (i\0) < CHEllill^llidlulli + KH0 



From (I4.26p . we now have, using the definition of negative norm 

(4.28) ||Et||_i < C||E||i + Ch\\p\\i + ^||E||i + Ch\\u t \\. 

Use Lemma 12.11 and Theorem 14.11 to complete the rest of the proof. □ 

Remark 4.1. Instead of using \J^.2J^ , if we have used a crude estimate, like ||Et||— < C||E t || (as 
in JIH (6.2)]), we have obtained 

In that sense, we have improved the pressure error estimate for conforming finite elements. 
Remark 4.2. In the process, we observe the following estimate, combining f^, 1\ ) and \J^.15 ): 



(4.29) e~ 2at f e 2as \\E(s)\\ 2 ds < Ke Kt h A . 

Jo 

Below we present an estimate for Et, which will be needed in our error analysis for two-level method. 
Theorem 4.3. Under the assumptions of Theorem \4-l\ we have 

(4.30) e~ 2at rai(s)||E s (s)|| 2 ds < Ke Kt h\ 

Jo 

Proof. The proof is similar to that of Theorem 14.11 and so, we will only sketch a proof here. First we 
split the error: 

E t = £t + »7f 

The equation in £ t is given by 

(4.31) (£ tt , cb h ) + uafa, <j> h ) = ( Pt , V • (f> h ). 

With P h : L 2 (f2) —> J h as L 2 -projection, we choose cb h = cri(t)P h ^ t = <Ji(t)(£ t — (u-P h u) t ) in (|4.31j) 
to find (recall a x (t) = (t* (t)) 2 e 2at , r*(t) = min{l,i}) 

^Mflll&ll 2 } + »*i(t)Mt\\l = l°u(t)Mtf + *i(f)«*, (u - p h u) t ) 

+vo-i(t)a(£ t , (u - P h u) t ) + ai{t)(p t - j h p t , V • Ph£ t )- 
12 



Use projection properties (Bl) and Cauchy-Schwarz inequality to obtain 

~{*i(«)ll* t H 2 } + "iWMtWl < Ca(t)Ut\\ 2 + ^|lK u " Wt\\ 2 
(4-32) +Ch.a 1 (t)\\£ t \\ 1 {\\u t \\ 2 + \\ Pt \\ 1 } 

(recall that a{t) = T*(t)e 2at , r*(t) = min{l,t}) 

Use Young's inequality and kickback argument. Integrate the resulting inequality to find 

<n(t)Mt\\ 2 + " f^WU^Wlds^c r*(8)\\Us)\\ 2 d8 + Ch.* 1 (t)\\u t \\l 

Jo Jo 

+Ch [ <T(s)\\u s ( S )\\lds + Ch 2 ! <Tl(s){\\u t \$ + \\ Pt \\l} d S . 

Jo Jo 

Use Lemma |2. II to observe that 

(4.33) {T*(t)) 2 \\£ t f + ve- 2at f ' a x {s)\\Us)\\lds < Ce~ 2at f a{s)\\£ s {s)\\ 2 ds + Kh 2 . 

Jo Jo 

To estimate the first term on the right hand-side of (|4.33p . we first recall the equation in 

(£t> 4>h) + <t>h) = (Pi v • 0h)> V< K e Jh- 
Choose cj> h = &(t)Ph£t above. 

*(t)Ut\\ 2 + ~{^(t)||€ll?} < Ce^uw 2 + cM*){ll€tlllMli + ll^llilKib + Mlill^lli}. 

Use kickback argument and then integrate with respect to time. Use Lemma 12. II to obtain 
Use (5.7) from [12], that is, 

(4.34) e~ 2at f e 2as U{s)\\\ds < Kh 2 

Jo 



to find 



(4.35) e~ 2at j a(s)\\Us)\\ 2 ds + r*(t)U\\l < Kh 2 + ^e~ 2at J a^U^Wlds. 
Using (Q51) in Bl3| leads us to 

(4.36) (r*(t)) 2 ||^|| 2 + e~ 2at f ai(s)\\Us)\\lds < Kh 2 . 

Jo 

Next, we use parabolic duality argument to establish (as in Lemma 5.1 from [T2], except that the 
right hand-side of (5.8), [12] should now read cri(t)£ t ) 

(4.37) e~ 2at f ai(s)U s (s)\\ 2 ds<Kh 4 . 

Jo 

An estimate of rj t would now complete the proof. By definition, we can easily deduce the equation 
satisfied by rj t . 

(4-38) (r] tt , <j> h ) + va{rj t , <f> h ) = A 1At ((j) h ), 
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where <t>\ht ls given by 

(4.39) Ai,ftt(0/J = -6(uh,t. E ' 0ft) " 6 K> E *' 0ft) - 6 ( E ^ u > 0ft) ~ 6 ( E , u <> 0ft)- 
Choose ft = (T\(t){— Ah)~ 1 r] t to find 

(4.40) |s{^(*)H*l*ll-i> + ^i(*)H»l*ll 2 ^ + oiWALfc.tCC-Afc)- 1 !!*). 
Using similar proof technique of (|4.14p . we estimate the non-linear term A-i h,t as follows: 

ffi(*)6(u fc)t , E, (-Afc) -1 ^) + <7i(f)6(E, u t , (— Aft) -1 T/ t ) 

<Cai(t)||E||||»| t ||i?||i| t || 1 / a (||u h , t ||i + ||ii t || 1 ) 

(4.41) ^(tJKI^ + Co-WH^IIii + Jfe^HEll 2 

and 

ffi(t)6(u fc> E t , (-Ah)" 1 ^) + <*\ (t)b(E t , u, (-Ah)" 1 ^) 
<Cfcl(t)(|&|| + II^IDhtll^ll^ll^dluftlli + Hulk) 
^i^ll^lllliltll^ll^ll^ + ^iWIIiitll^lliitll 8 ^ 

(4.42) ^WII^f + ^WII^H^ + KaxWII^II 2 . 

Incorporate (I4.41|) and (|4.42p in (|4.40p . Use Young's inequality and appropriate e > and finally, 
integrate with respect to time to find 

vi(t)\\Vt(t)\\-l + >> f <Ti(s)\\r, 8 (s)fds < K [ t <T(8)\\r la (8)\\ 2 _ 1 ds + K f e 2as \\E(s)fds 
Jo Jo Jo 

+K [ a^WUsWds. 
Jo 

Use (|4T29l) and (Q71) to obtain 



t rt 

Ktu* 



(4.43) *i(*)Ik(*)||ii + W a 1 {s)\\r 1s {s)fds<K{t)e Kt h^K a{s)\\n s {s)\\i x ds. 

Jo Jo 

To estimate the last term of (|4.43p . we recall the equation in rj (|4.t>[) : 

(Vt> 0ft) + ua (V, 0ft) = A i,ft(0ft)- 
Choose 4> h = a(t)(—Af l )~ 1 rj t to find 

(4.44) *(*)ll»Itl|li + ~{^(*)||»»|| a } < ^t(t)\\vf + .(^^((-Aft)- 1 ^. 

Similar to (|4.14p . we estimate the non-linear term as follows: 

a(i)Ai )/l ((-A ft )- 1 » 7t ) = -a(t)b(B, u h , (-A^Vt) ~ <r(t)b(u, E, (-A*)" 1 ^) 
<C^)||E||h t ||_ 1 {||u,|| 2 + ||u|| 2 } 

Incorporate this in (|4.44p . use kickback argument and then integrate with respect to time. Next, use 
Lemmas 12.11 and 13.21 

(4.45) [* v(s)\\vM\\-lds + cx(t)\\v(t)\\ 2 <C [' e 2as \\v(s)fds + K f e 2as \\E(s)fds. 
Jo Jo Jo 
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Apply (|4.15p and (|4.29|) in (|4.45p . Use the resulting estimate in (|4.43j) to obtain 
(4.46) ^WINtWf-i + * /"ffi(s)|to»[| a cfe < K(t)e Kt h\ 



This along with ()4,37p now completes the rest of the proof. □ 
The pressure error estimate for the two-level method requires another estimate. 
Lemma 4.2. Under the assumptions of Theorem \4.1[ we have 

(4.47) 



! a 1 {s)\\E s {s)f 1 ds<Ke Ki h 2 . 
Jo 



Proof. Keeping in mind (14.36[) . we only need the estimate in rj. Choosing 4> h = ai(t)r) t in (I4.39P and 
estimating non-linear as 

A w fa t ) < s\\ Vt \\l + C||E||f (\\u htt \\i + HHx) + C||£ t ||?(H||? + ||u||?) + ^II^H 2 , 

we find 



<ri{t)\\Vt(t)\f + " *i(s)\\ri.(s)\\ids<K o-(s)\\ Vs (s)\\ 2 ds + K a^UMWtds 
Jo Jo Jo 

+Ke Kt h 2 ra( S ){||u M ( S )||? + ||u s ( S )||?}d S 
Jo 

Using (|4.36p and Lemmas 12.11 and 13.31 we have 

(4.48) a 1 (t)\\ Vt (t)f + ^ ! <Ti{s)\\r)s{s)\\lds<K f cr( S )\\r, s (s)\\ 2 ds + Ke Kt h 2 . 

Jo Jo 



Now, for the first term on the right-hand side, we take <p h = cr(t)rj t in the rj equation, to obtain, as 
in (Q4"j) : 

°(t)\\Vt\\ 2 + !^M*)NI?} < ^t(t)\\v\\l + *(t)A 1>h (rj t ). 
We estimate the non-linear term as: 

A 1) ,(T 7t )<^||r 7t ||||E|| 1 (||u /l || 2 1/2 + ||u||f ) 

to find 



(4.49) 



°(t)\\v(t)\\l + f a{s)\\ Vs {s)\\ 2 ds<Ke Kt h 2 + C f e 2 ^\\r,( S )\\ 2 ds. 
Jo Jo 



To estimate the last term, choose <f> h = e 2at r] in the rj equation and estimate the non-linear term as: 

Ai,/ l (r 7 )< e ||r 7 !| 1 +A'(^|| 2 + h|| 2 ) 

to obtain 

|^l| 2 + ^lli<^ll? + ll^l| 2 ) 
Integrate with respect to time and the use (5.7), [12J, that is 

(4.50) I' e 2as U{s)\\jds < Ke 2at h 2 , 

Jo 
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to observe 



\f)(t)f+ I' e 2as \\v(s)\\jds < Ke 2at h 2 + I e 2as \\r]{S) \\ 2 ds. 
Jo Jo 



Use Gronwall's lemma. 



(4.51) II^WH 2 + I' e 2as \\v{s)\\ids < Ke Kt h 2 . 

Jo 

Combining (|4.50p and (|4.5ip . we conclude 

„2asi|-nV„M|2 . ^ T s„Ktr2 



(4.52) / e zas \\B{s)\\ids < Ke^h 1 

Now from (|4.49j) . we have 



o 



(4.53) v(t)\\v(t)\\l+ f 'a(s)\\r) s (s)\\ 2 ds<Ke Kt h 2 . 



o 



And from (|4.48|) . we have 



(4.54) <n(t)\\Vt(t)f + " f o^W^Wlds < Ke Kt h 2 . 

Jo 

Combining with (|4,36p , we complete the rest of the proof. □ 

5 Two-Level FE Galerkin Method 

In this section, we work with an additional space discretizing parameter H, that corresponds to a 
coarse mesh. In other words, < h < H and both h and H tend to 0. We introduce associated 
conforming finite element spaces (Hh,Lh) and (H/j,L^) such that (H#,L#) C (H^L^). And this 
two- level finite element is to find a pair (u h ,p h ) as follows: 

First Level: We compute the mixed finite element approximation (uh,Ph) £ (H-h,Ph) of (u,p) 
of (I2.2p . In other words, we solve the nonlinear problem in a coarse mesh. Find (uh,ph) € (Hh,Ph) 
satisfying 

(iifft, 4> H ) + va(u H , <f> H )+b(u H ,u H , <f> H ) - (p H , V • <f> H ) = (f , <p H ), 

(5.1) (V-u h , X h) = 0, 

for {4>h,Xh) S {il H ,L H ). 

Second Level: We solve a linearized problem on a fine mesh. In other words, we solve a Stokes 
problem. Find (u h ,p h ) £ (H^L^) satisfying 

("t > <t>h) + va(u\ <t>h) ~ (P H , V-<f>h) = ( f > <t>h) - K^H, u H , <f> h ), 

(5.2) (V-u h , X h) = 0, 

for {4>hiXh) 6 (H h ,L/j). 

An equivalent way is to look for solution in a weekly divergent free space. 
First Level: Find u# € J// satisfying 

(5.3) (u m , 4> H ) + ra(u fl , 0# ) + b(u H , u H , <j> H ) = (f , <f> H ), 
for c() H eJ H . 

Second Level: With u# as the solution of (15.31) . find u ft 6 J/j satisfying 

(5.4) (u t h , 4> h ) + ua{u h , 4> h ) = (f, 4> h ) - b(u H , u H , 
for <t> h €J h . 
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Remark 5.1. The well-posedness of the above systems can be seen from the facts that \5. 1\) or \5.3\) 
is classical Galerkin approximation and hence is well-posed as is stated in Section 3. And \5. Sty or 
( 5.4\ ) represent linearized version. Given ujj and with suitable u h (0), it is therefore an easy task to 
show the existence of an unique solution pair (u h ,p h ) (or an unique solution u h ) for the linearized 
problem following the foot-steps of the non-linear problem. 

Following Lemma 13.21 we can easily obtain the a priori estimates of u#. 

Lemma 5.1. Under the assumptions of Lemma \3.S\ and for u;/(0) = PhUq, the solution u# of \5.3\) 
satisfies, for t > 0, 

(5.5) \\u H (t)\\ +e~ 2at f e 2as \\u H (t)\\ 2 ds < K, 

Jo 

(5.6) ||iitf(0||i + e- 2a * f e 2as \\u H (t)\\ 2 ds < K, 

Jo 

(5-7) (T*(t)y/ 2 \\u H (t)\\ 2 <K, 

where r*(t) = min{l,t} and K depends only on the given data. In particular, K is independent of 
H and t. 

The following higher-order estimate of u# is required for error analysis. The proof of the same is 
similar to that of Lemma [ 



Lemma 5.2. Under the assumptions of Lemma \5.1\ the solution u# of $5.3}) satisfies, fort > 0, 

(5.8) (T*) 2 (t)\\u H , t (t)\\l + (r*f(t)\\u H , t (t)\\ 2 + e~ 2at f T*(s)e 2 ™\\u H>s (s)f ds < K. 

Jo 

For the a priori estimates of u h , we present another lemma. 

Lemma 5.3. Under the assumptions of Lemma \5.1\ the solution u h of \5.J^ satisfies, for t > 0, 

(5.9) \\u h (t)\\+e- 2at [* e 2as \\u h (t)\\l ds < K, 

Jo 

(5.10) ||u h (i)||i + e~ 2at [ e 2as \\u h (t)\\ 2 ds < K, 

Jo 

Proof. Given u# along with the estimates of Lemma (|5.ip . we choose 4>h = e 2at u h (t) = e at u h (t) in 
(foTil) to obtain 

(5-11) ^II^H 2 " "H^ll 2 + "11**11? = (f , u h ) ~ e 2at b(u H , u H ,u h ). 

Use Cauchy-Schwarz inequality, Poincare inequality with first eigenvalue of Stokes operator as the 
constant and Young's inequality, we have 

(5-12) &u h ) < [|f||[|fi h || < n U|f||||u' l || 1 < Un h \\\ + -^-\\f\\ 2 . 

A} 7 4 v\i 

From Lemmas 13.11 and 15.11 we find that 

6(u / ,,UH,^)<||u / ,|| 1 /2|| UH ||f || u& ||l/2|| u& ||l/2 + || u ^|||| UH || 1 ||u fc || 1 

<K+ U -\\A\l 
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Therefore, 

(5.13) e 2at b(u H ,u H ,u. h ) < Ke 2at + ^\\u h \\ 2 v 



Putting the estimates f|5. 12j) - (|5. 13 j) in (|5.1ip gives us 



^\\u h \\ 2 - 2a\\u h \\ 2 + u\\u h \\ 2 < Ke 2at + -i-||f II 2 . 
Use Poincare inequality to get 

(5.14) ±\\u h \\ 2 + \\u h \\ 2 < Ke 2at + -U|f || 2 . 

at Ai i^Ai 

Since < a < uXi/2, we have that z/ — 2a/ Xi > 0. 
Integrate (|5.14p with respect to time. 

\\u h (t)\\ 2 + („ - ^) T [|ft h ( a )||?ds < Ke 2at . 
Ai Jo 

Multiply by e" 2crf to conclude (pT9l) . 

For the next estimate, we choose cf) h = e 2at AhU h (t) in (|5.4|) and proceed as above. For the non-linear 
term, we again follow similar proof technique of (|4.14p to obtain 

(5.15) b(u H ,u H ,A h u h ) < C||u // || 1 / 2 ||u jF/ ||2 /2 ||u//||il|u 
These estimates let us have 



2- 



(5.16) i^Hu^ii? — ciiu^nf -k < ||F|| ||u^|| 2 -k Ce^llu^ll 1 / 2 !!^^!^/ 2 !!^^!!!!!^!^. 

Using Cauchy-Schwarz inequality and Poincare inequality as earlier, we find that 

(5-17) ^\\u h \\l + {» ~ ^) ll^lll < C\\f || 2 + Ce 2at {\\u H \\ 2 \\u H \\i + ||u„|| 2 }. 

Integrate (I5.17P with respect to time, use Lemma ED and finally multiply by e~ 2at to obtain (I5.10p . 



6 Error Estimate 

In this section, we present the error estimate for the spatial approximation, that is, two-level finite 

element approximation. We achieve the desired results through a series of Lemmas. For the sake of 

convenience, henceforward, we will write Ke Kt simply as K. 

We denote the error, due to two- level method, as e = — u h . 

From the equations (13. 31) and (15. 4h . we have the following error equation: 

(6.1) {e t ,4> h ) + va(e,4) h ) = K h (4> h ) V<j> h G 3 h , 
where 

^h{4>h) = b ( u H, Uff, 4>h) ~ H u h, u h , <t> h ) 

(6.2) = b(u H , u H - u h , 4> h ) + b(u H - lift, u h , 4> h ). 
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Lemma 6.1. Under the assumptions of Theorem \^.l\ and with the additional assumption that u 
11^(0) and for < a < u\\/2, the following estimate 

(6.3) \\e{t)\W + e- 2at f e 2as \\e{s)\\ 2 ds < K H A 

Jo 

holds, for t > 0. 

Proof. Choose 4> h = e 2at {- K h )- 1 e(t) = e at {- K h )- 1 e{t) in (HH) to obtain 

(6-4) \jt M2 - 1 ~ aM - 1 + = ^^({-^hT 1 *). 

Following (|4.14p and using Lemmas 13.21 and 15.11 we find 

A((-A h )- 1 e) = b(u H , u H - u h , (-A h ) _1 e) + b{u H - u h , u h , (-A/ l ) _1 e) 
< ||uff - u^||||e||{||uH||i + \\uh\\i} < K\\u H - u^||||e||. 

Incorporate this in (|6.4p . 

(6-5) Jt^- 1 + ~ 2a ^- 1 + K ^ H ~ ilh ^- 

Integrate (|6.5p with respect to time. Using (|4.29p . we first observe that 
«t rt r-t 

(6.6) 

And hence 



[ \\(u H -u h )(s)\\ 2 ds< [ \\(u-u H )( s )\\ 2 ds+ [ \\(u-u h )(s)\\ 2 d S <Ke 2at H A . 
Jo Jo Jo 



(6.7) \\e(t)\\-i + » f l|e(s)|| 2 ds < 2a f Ms^ds + Ke 2at H A . 

Jo Jo 



'0 Jo 

-2at 



Multiply by e and use Gronwall's Lemma to complete the rest of the proof. 
Lemma 6.2. Under the assumptions of Lemma \6.1\ the following estimate 

(6.8) r*(t)||e(t)|| 2 + e~ 2at f a(s)\\e(s)\\ 2 ds < K(t)H 4 

Jo 

holds, for t > 0, where a(t) = T*(t)e 2at , r*(t) = min{l,t}. 
Proof. Choose <p h = a(t)e(t) in (|6.ip to obtain 

(6.9) ±{ a ( t )\\e\\ 2 } + 2v a(i)||e|| 2 = ^(t)||e|| 2 + 2a(t)A h (e). 



To estimate the non-linear term, we follow (|4.19|) and use Lemmas 15.11 and 13.21 

Afc(e) < C\\u H - ||e||i{ ||u jF/ ||i /2 ||ujy||2 /2 + ||u/ l ||i /2 ||u/ l ||2 /2 } 
(6-10) < K\\u H - u fc ||||e||i{||u H ||2 /2 + \\u h \\l /2 }. 

And hence, we find from (16.91) 



^{a(t)\\e\\ 2 } + 2v a(t)\\e\\ 2 < C\\e\\ 2 + ^(t)||e||i||u H - u h \\{\\u H \\l /2 + K||* /2 }. 
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Use Cauchy-Schwarz inequality and then kickback argument to yield 

(6.11) ±{a(t)\\e\\ 2 } + v a(t)||e||2 < C||e|| 2 + Ka(t)\\u H - u h \\ 2 { \\u H \\ 2 + \\u h \\ 2 }. 
Use Lemmas 15.11 and 13.21 to estimate the last term of (|6,lip as follows: 

Kcr(t)\\vL H - u h \\ 2 {\\\i H \\ 2 + ||u ft || 2 } < K\\u H - u h \\ 2 . 

We have used the fact that T*(t) < 1. Incorporate this in (|6.1ip and integrate in time. Use (|6.6p to 
estimate the last term. 

(6.12) a(t)\\e(t)\\ 2 + v f a(s)\\e(s)\\lds <C I \\e{s)\\ 2 ds + Ke 2at H\ 

Jo Jo 

Now use Lemma 16. II and then multiply the resulting inequality by e~ 2at to complete the rest of the 
proof. □ 

Lemma 6.3. Under the assumptions of Lemma \6.1\ the following estimate 

(6.13) (T*(i)) 2 ||e(i)||? + e- 2Qt f ^(s^e^fds < K(t)H 4 

Jo 

holds, fort>0, where a x {t) = (t* (t)) 2 e 2at . 
Proof. Choose <f> h = oi{t)e t {t) in (|6.ip to obtain 

(6.14) 2a 1 (t)||e t || 2 + i/^{<7i(*)||e||?} = ^ ljt (t)||e|| 2 + 2a 1 {t)K h {e t ). 
Integrate (|6.14p with respect to time. 

(6.15) 2/ ai(s)||e s (s)|| 2 (is + ^i(t)||e(t)|| 2 < C I a(s)\\e(s)\\ 2 ds + 2 f ax(s)A h (e s ) ds. 

Jo Jo Jo 

We rewrite the nondinear terms as follows: 

A h (e t ) = b(u H , u H - u h ,e t ) + b(u H - u h , u h , e t ) 

= ^{b(u H ,u H - u h ,e) + b(u H - u^u^e)} - b(u H ,t,u H - u h ,e) 
- b(u H , u H ,t ~ u h ,t, e) - b(u H ,t - u hit , u h , e) - b(u H - u h , u h>t , e). 

And hence 

eri(i)Afc(e t ) = cri(i)(&(uff,iiH - u h ,e t ) + b(u H - u h ,u h ,e t )) 

= -^{o-^^uh^h - u h ,e) + b(u H - u h ,u h ,e))} 

- a ljt {t)(b(u H ,u H - u h ,e) + b(u H - u h ,u h ,e)) - ax(t)b(M H ,t, u H - u h ,e) 

(6.16) - ai(t)(b(u H ,u H ,t ~ u htt ,e) + b(xi H)t - u htt ,u h ,e) + b(u H - u h ,u h<t ,e)). 

As seen earlier, we have 

-b(u H ,t, u H - u h ,e) - b(u H - u h , u h;t ,e) 

, ^ u r 1 1 11 1/2 n ,, 1/2 1/2 ,, ..1/2 1 

< C\\u H - Uh||||e||i{||uH,t||i ||u^|| 2 ' +11^^11/ Hu^illa' ) 
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Use Lemmas 13.31 and 15.21 to conclude that 

(6.17) - b(u Hit ,u H - u h ,e) - b(u H - u h ,u hjt ,e) < K(T*(t)y 3/2 \\u H - u h \\\\e\\i. 
Similarly 

(6.18) - b(u H ,u H - u h ,e) - b(u H - u h ,u h ,e) < K(t* (t)y 1/2 \\u H - u h \\\\e\\x 
and 

(6.19) - b(u H ,u Hit - u M ,e) - b(u H j - u h j,u h ,e) < K(r*(t))~ 1/2 \\u Hjt - u M ||||e||i. 
Incorporate (|6.17p - (|6.19p in (|6.16p and integrate with respect to time. Re-use (16.18P to find 

/ ai(s)A h (e s ) ds < ai(t)(b(u H ,u H - u h ,e) + b(u H - u h ,u h ,e)) 
Jo 

Jo 

<Ke 2at (r*(t)f/ 2 \\u H -u h \\\\e\\ 1 + C f a(s)\\e(s)\\ 2 ds 

Jo 

+ K fe^iUuH-Uh^W 2 + (r*( S )) 2 ||(u^-u v )( S )|| 2 } ds 
Jo 

^amnl + KamuH-Uhf + C f a(s)\\e(s)\\ 2 ds 
1 Jo 

(6-20) +K r e 2 -{||(uH-u h )( S )|| 2 + (r*( S )) 2 ||(u^-u M )( S )|| 2 } ds. 

Jo 

Put (f6T20]) in ([67T5]) to obtain 

2 f cn^lle,^)!! 2 ^ + £ai(t)||e||? < C f a( S )\\e(s)\\ 2 ds + Ka(t)\\u H - u h \\ 2 
Jo 1 Jo 

(6-21) +K re^jlKu^-uO^f + ^^^IKu^-u^)^)!! 2 } ds. 

Jo 

Now, use Lemma 16.21 For the remaining part, rewrite u# - uj, as u - u/, - (u - u#). Then use 
(I4.29P and Theorems 14.11 and 14.31 to complete the rest of the proof. □ 

Lemma 6.4. Under the assumptions of Lemma \ 6.1\ and additionally that the assumption (B2') 
holds, we have 

(6.22) \\Ph-p h \\ < i^XT^))- 1 ^ 2 . 

Proof. The LBB condition (B2') tells us that, for t > 

hn _ (Ph-P H ,V -ct> h ) 



(6-23) \\Ph-F\\i?f*<Ko SU P 



<#0 h eH h II Vh 111 



From (|3.2p and (|5.2p . we have, for <fi h £ 



(Ph - P h , V • 4> h ) = {e t , cf> h ) + ua(e, <j> h ) - A h ((f> h ) 

< Cll^lhllletH^ + Helix} + K(T*{t))- x '*\\viH ~ Vh\\\\<t> h \\i- 
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We have estimated as in (|6.10p and have used Lemmas 15.11 and 15.31 Using (|6.23p . we obtain, 
thanks to Lemma 14.11 and 16.31 



(6-24) \\Ph-p h \\ L yR <C\\e t \\^ lth + K{T*{t))- l H 2 . 

Following (|4.24j) . we observe 

||et||-i,/i < ||et||_i. 

Keeping in mind that, for (f> £ Hq, Ph4> G J/u we use (|6.ip to write (e^, <f>) as in (|4.26p and similar to 
estimates in (|4.27p . we find 

(6.25) 

(et, 0) < C||e||i||0||i + C\\u H ~ u h \\{\\u H \\{ /2 \\u H \\ l 2 /2 + K||; /2 K||J /2 }|M|i + Ch\\e t \\ ||0||i. 

Use Lemmas 15.11 15.31 and 16.31 We again rewrite uh — as u — — (u — u#) and use Theorem 14.11 
Now, for ^ 0, we divide (I6.25P by to finally obtain 

HetH^^r*^))- 1 /? 2 + Ch\\et\\. 

And hence, from (|6.24p . we find 

\\Ph " All?/* < CH\\e t \\ + K^it))- 1 !! 2 . 
Now the following lemma completes the proof. □ 
Lemma 6.5. Under the assumptions of Lemma \6.1\ the following estimate holds, for t > 0. 

IM| <K{r*{ft)- x H. 
Proof. Differentiate the error equation (|6.ip with respect to time. 
(6-26) (e tt ,(f> h ) + i'a(e t ,(f> h )=A hit ((f) h ) V<f> h € 3 h , 

where 

^h,t{<t>h) =K u H,t, uh - u h , 4> h ) + b{u H , u H ,t ~ u ft)t , 4> h ) + b(u H ,t ~ vi hit , u h , 4> h ) 

(6.27) + b(u H - u h ,u htt ,(f) h ). 

Choose cj> h = a\(t)e t in (|6.26[) to find 

(6.28) ^ W*)l|e 4 || 2 } + 2i/<7 1 (t)||e t || 2 < a M (t)||e t || 2 + 2a 1 {t)A h>t (e t ). 
Similar to (|6.17p - (|6.19p . we obtain 

2<r 1 (t)A htt (e t ) < Cax(t)\\e t \\i{\\u H ,t\\i + \\u h j\\i}\\u H - u h \\i 

(6.29) +Ccri(t)||et||i{||u/f||i + HuhHijUu^t - u h>t \\i. 

Incorporate (|6.29p in (|6.28p . Use kickback argument and then integrate with respect to time. 

a 1 (t)||e 4 (t)|| 2 + ^ f* <Ti(s)\\et\\l <C f a(s)\\e t (s)\\ 2 ds + K I ai{s)\\{u HjS - u M )(a)||?da 
Jo Jo Jo 

(6-30) +KH 2 [ a(s)(\\uH, s (s)\\l + K, S ( S )H?) 

Jo 
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Use Lemmas 12.11 and 15.31 to bound the last term of (|6.30p . For the second last term, we rewrite 
uh )S — ^h,s as u s — — (u s — uh, s ) and use Lemma [4^21 



(6.31) <7 1 (i)||e t (i)f + I , f ai(s)\\e t \\i < C f a(s)\\e t {s)\\ 2 ds + KH 2 . 

Jo Jo 

To estimate the first term on the right hand-side of (|6.31|) . we put 4> h = cr(t)e t in 

(6-32) ^Wlletll 2 + ^M*)H e ll?} = <rt(t)\\*\\i + *(*)A h (et). 

We take care of the non-linear term as earlier. 

A h (e t ) = b(u H , u H - u h , e t ) + b(u H - u h , u h , e t ) 

< C\\e t \\\\u H - UfcllidluffH^HuHllJ/ 2 + KII^KI^ 2 ) 

+ C||e t ||||uH - iVj|| 1/2 ||uh - u^|} /2 (||u#||i /2 ||uh||2 /2 + Hu/, || } /2 ||u ft H^ 2 ) 

< K\\e t \\\\u H - u h ||i([|u ff ||^ /2 + K||^ /2 ). 

We have used Lemmas 15.11 and I3.2L Incorporate in (|6.32j) and after kickback argument, integrate 
with respect to time. 

t ft i-t 

2j„ , _AAII„AAIl2 ^ I „2as\\„f„\\\2 j„ i v / „2asi|T7i/„M|2, 



a(s)\\e s (sWds + a(t)\\e(t)\\i<C / e ||e(s)||fds + K / e 2as \\E{s)\\{ds. 
o Jo Jo 

Use (g32D to obtain 

(6.33) f a( S )\\e s (s)fds + a(t)\\e(t)\\ 2 <C f e 2as \\e{ S )\\\ds + KH 2 . 
Jo Jo 

To estimate the first term on the right hand-side of (|6.33p . we choose 4> h = e 2at e in (|6.ip to find 

\j t \\e\\ 2 + [v - j-) \\e\\ 2 < e 2at A h (e) < ^^^{e^Hua - u ft ||i}. 

Use kickback argument and then integrate with respect to time. 

e 2at ||e(t)|| 2 + f e 2as \\e{s)\\\ds<K f e 2as \\E{s)\\\ds. 
Jo Jo 

We recall that e(0) = 0. Apply P~52"]) to get 

(6.34) e 2at \\e(t)\\ 2 + f e 2as \\e{s)\\lds < KH 2 . 

Jo 

Now, from (|6.33p . we have 

(6.35) I a{s)\\e s {s)\\ 2 ds + a(t)\\e{t)\\j < KH 2 . 

Jo 

And from (gSjTJ 

(6.36) ffi(*)||et(t)|| 2 + i/ r^Wlletll? <KH 2 . 

Jo 

And this completes the rest of the proof. □ 
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Remark 6.1. In the Lemma \6.5l we have desired and obtained a lower order estimate, which is of 
0(H) instead of 0(H 2 ). Note that in the proof, we have avoided 0(H 2 ) estimates. For example, for 
the second last term &6. °JU\) . we do not follow the argument that we have used to bound the last term 
of i6.21\) . Instead, we have used Lemma \4-S\ An optimal estimate of Ije^l would read 0(t~ 3 ^ 2 H 2 ) 
and which could be achieved by choosing 4> h = o~2(t) e t i> n lib' -26}) . instead of o~i(t)e t as in the proof 
of Lemma \6.b\ But this would lead to a higher order singularity of pressure error estimate at t = 0. 
Our pressure error estimate \6. 22\) and velocity error estimate \6.13\) have same order of singularity 
att = 0, which is not the case in f^jj. Even for smooth initial data, we can obtain error estimates for 
both velocity and pressure, that have same order of singularity at t = 0, i.e, of 0(H 2 ), by following 
our proof technique. This improves the pressure error estimate of J2[/. To be precise, in J2J/ ; to 
estimate the pressure error, ||ej|| is used (see (5.39) from JE/). Taking advantage of conforming finite 
element, we have used 1 1 e £ | { i instead. And we are awarded with a better pressure error estimate. 
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